Double non-equivalent chain structure on vicinal Si(557)-Au surface 
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We study electronic and topographic properties of the vicinal Si(557)-Au surface using scanning 
tunneling microscopy and reflection of high energy electron diffraction technique. STM data reveal 
double wire structures along terraces. Moreover behavior of the voltage dependent STM tip - surface 
distance is different in different chains. While the one chain shows oscillations of the distance which 
are sensitive to the sign of the voltage bias, the oscillations in the other chain remain unchanged 
with respect to the positive/negative biases. This suggests that one wire has metallic character 
while the other one - semiconducting. The experimental results are supplemented by theoretical 
calculations within tight binding model suggesting that the observed chains are made of different 
materials, one is gold and the other one is silicon chain. 
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I. INTRODUCTION 

Studies of the high-index (vicinal) surfaces have at- 
tracted much attention as templates for formation of the 
low dimensional structures on them Q. Such structures, 
usually one-dimensional, are important from scientific 
point of view as they allow to study phenomena like Lut- 
tinger liquid H0, Peierls instability Q, self-assembling 
of arrays of nanowires 5] or basic properties of electrons. 
On the other hand they have potential technological ap- 
plications in nanoelectronics or quantum computing. 

There are many examples of the vicinal surfaces stud- 
ied recently, including Au decorated Si(335) surface 
clean Si[557) 0, Si(557)-Au H El El 
111 HQ, Si(557)-Ag or Si(5512)-Au 
l23| . All these surfaces allow to grow on them 
one dimensional wires along terraces. 

In the case of the Si(335)-Au surface, which con- 
sists of Si(lll) terraces 3| x wide, a single long 
mono-atomic chain on each terrace has been observed 
at gold coverage 0.28 ML. This surface has been inves- 
tigated using reflection of high energy electron diffrac- 
tion (RHEED) angle resolved photoemission spec- 
troscopy (ARPES) (3, 0| and scanning tunneling mi- 
croscopy (STM) indicating metallic character of the 
wires grown on it. Moreover, it was found experimentally 
and confirmed theoretically that these wires show bias 
dependent oscillations of the STM tip - surface distance 
along wire with a period of 2 x a^jo] [i| 

Jalochowski et al. showed how small amount of 
gold stabilizes Si(557)-Au surface. This surface possesses 
wider terraces (5| x a^j) and at gold coverage 0.2 ML 
one observes not single but double mono-atomic chains 
on each terrace 

111 El El 111 El El El El Many ex- 
perimental techniques have been applied to study prop- 
erties of this surface. Segovia et al. |l2| using ARPES 



claimed that the wires show Luttinger liquid (LL) be- 
havior which turned out later on to be not LL but usual 
Fermi liquid with two metallic bands [IH . . Those con- 
clusions have been reached by studying it using ARPES 
and STM techniques and confirmed theoretically within 
density functional theory 15]. On the other hand, con- 
ductivity measurements suggest non-metallic character 
of the surface which seems to be in contradiction 
with ARPES data. However, one has to have in mind 
that the chains are not ideal but they are discontinued 
(tens of nanometers long) so even along wire there are 
metallic as well as semiconducting regions, which may 
explain discrepancy between conductivity measurement 
and ARPES data. Moreover, it turned out that such 
amount of gold is sufficient to form single wire only 
H ^3 •> unless we assume that gold atoms are distributed 
with distance equal to twice of the lattice constant in 
direction [110]. Thus it is enough of gold to form two 
chains. However the existence of the single Au chain per 
terrace has been found using x-ray diffraction p^ . 

In the present work we use scanning tunneling mi- 
croscopy (STM) technique to study gold induced struc- 
tures on Si(557) surface. Those experimental studies are 
supplemented by theoretical ones based on tight binding 
model where the STM tunneling current is calculated 
using non-equilibrium Keldysh Green function formal- 
ism. Theoretical description remains in good qualitative 
agreement with the STM data. Rest of the paper is orga- 
nized as follows: in Sec. [H]we describe the experimental 
setup and provide some data. In Sec. IIIII we introduce 
model of the chains on surface and the results are pre- 
sented in Sec. IIVI The comparison of the theoretical 
calculations with the experimental data is given in Sec. 
IVl Finally Sec. IVII contains some conclusions. 



II. 



EXPERIMENTAL 
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Experimental setup consists of ultra high vacuum 
(UHV) chamber equipped with a scanning tunneling 
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microscope (type OmicronVT) and reflection high en- 
ergy electron diffraction (RHEED) apparatus. Samples 
were prepared in-situ and the base pressure was less 
than 5 x 10 -11 mbar during measurements. The one- 
dimensional structures have grown after deposition of 0.2 
ML of Au, heating the sample at temperature 950 K for 
20 s and gradually annealing to the room temperature 
for 3 min. The quality of the surface reconstruction has 
been controlled by RHEED technique. All STM data 
have been collected at room temperature. 

Figure E] shows large area STM image of the Si(557) 
surface covered with 0.2 ML of gold. 




FIG. 1: (Color online) The 30 x 30 nm 2 area STM image 
of the Si(557) - 0.2 ML Au surface recorded with tunneling 
current I = 0.05 nA and sample bias equal to U = — 1.0 V. 
The frames mark areas chosen for further analysis. 

Figure shows enlarged part of the image marked with 
thick line in Fig. n recor ded at tunneling current / = 0.05 
nA and the sample bias U = —1.0 V (a) and +1.0 V (b). 
We note that images (a) and (b) represent the same area 
of the sample. The STM images clearly show sets of dou- 
ble chains on each terrace. Distance between these sets is 
equal to the width of the unreconstructed Si(557) surface 
terrace, 5| x a^2] (1.92 nm), and the distance between 
two chains on the same terrace is equal to 2 x a^n2] (0-68 
nm). However the chains on the same terrace show dif- 
ferent behavior. While the chain indicated by arrows and 
denoted as C in Fig0 (a) and (b) is apparently continu- 
ous and very similar to that grown on Si(335)-Au surface 
the wire D differs considerably and reminds a row 
of well ordered single atoms rather than one-dimensional 
object. However periodicity along chains are the same in 
both wires and equal to 2 x a^jo] (0.77 nm). It is worth 
to note that similar periodicity along the chain has been 
also found in different structures, like Si(335)-Au 9] with 
single chain per terrace or Si(lll)-(5 x 2)-Au 24] with 
double Au chains. 

Figure |3| displays cross sections along atomic chain C 




FIG. 2: (Color online) The 9.75 x 13.5 nm 2 STM topography 
images of the same area of the Si(557)-Au surface recorded at 
two different sample biases U = — 1.0 V (a) and U = 1.0 V 
(b) with the tunneling current / = 0.05 nA. Bottom panels 
show profile lines perpendicular to the chains, indicated by 
long arrows in the main panels. 



(a) and D (b) indicated by arrows in FigEl As one can 
read from the figure both chains show oscillations of the 
electron density with distance along their length. More- 
over, it is clearly seen that those chains behave differ- 
ently upon reversing the sample bias. Topography of the 
nanowire C (top picture) reverses (maxima change into 
minima and vice versa) with switching of the voltage bias 
from negative to positive values. Such behavior has been 
previously observed in the case of Si(335)-Au structure 
|£j and has strictly electronic origin. On the other hand 
there is a lack of such reversing in the chain D where its 
topography remains almost unchanged, except that at 
positive bias the structures are higher. We stress that in 
order to exclude possible influence of the sample thermal 
drift the presented STM images were recorded simultane- 
ously, with negatively biased sample during line forward 
scan and positively biased during line backward scan. 

It is interesting to see the changes of topography with 
the bias voltage. The evolution of topography changes 
with the bias voltage is shown in Fig. 21 We observe 
that the chains of C (D) type simultaneously become 
less (more) visible while going from negative to positive 
biases. Moreover, chains of D type do not show any to- 
pography reverse. On the other hand, there is a reverse 
of the topography of the C type chains when crossing 
zero bias. 

Another interesting effect is associated with recently 
observed zero dimensional end states 25]. Such states 
must have the wave function localized to the end atoms 
in a chain, and it must must decay exponentially into 
the chain IH Similarly for Si(553)-Au surface 25], 
we observe that C type chains at negative sample bias 
appear longer than at positive bias (see Fig. |3| - top 
panel). At negative bias the end atoms are well visible 
while for negative bias the atoms second form the end 



3 



E 







; V \/ \i\f \i\ !\ / wVH l'/ 

chain C 








U = -1.0 V - 








1.0 V 





2 4 6 8 10 12 



nm 

0.22 | ■ ■ ■ 

j, 

0.20 •!• Hi |S i A A A , 




2 4 6 8 10 12 
nm 



FIG. 3: (Color online) Cross sections along middle chains 
indicated by short arrows in Fig. 0(a) and (b). Top panel 
shows cross section along chain C for bias voltage U = — 1.0 
V (solid line) and +1.0 V (dashed line), the bottom panel - 
along chain D. 



are enhanced and the end atoms became supressed. This 
fact can be explained in terms of different densities of 
states above and below the Fermi energy (see discussion 
in Sec. HVj) . However, unlike in Ref. |25j . the topography 
modulation along chain is twice a lattice constant (2 x 
a [iIo] =0.77 nm) and what is important we also observe 
the topography reverse. We stress that there is no such 
an effect for chains of D type (see Fig. bottom panel). 

To see the evolution of the topography of chain C type 
with the sample bias we show in Fig. 0the cross sections 
along the longest chain indicated by arrow with label C 
in Fig. [3 It is clearly visible how the length of the chain 
becomes shorter going from negative to positive voltages. 
Moreover around zero bias we observe the reverse of the 
topography. 

Figure shows the length of above discussed the chain 
as a function of the sample bias. We define the length as 
a distance between endmost maxima in the topography 
profile (see Fig. 0). We observe that the effect associated 
with the length changing of the chain takes place in the 
range between —0.75 V to +0.5 V. Below the bias —0.75 
V and above +0.5 V it remains almost constant, and the 
length difference between those biases is equal 0.77 nm 
(twice a lattice constant along the chain). 

From the experimental data we cannot judge what ma- 
terial (Si or Au) those wires are composed of. Certainly 
there is no sufficient amount of gold to produce two gold 
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FIG. 4: (Color online) The 9.6 x 19.1 nm 2 STM topography 
images of the same area of the Si(557)-Au surface recorded at 
different sample biases (indicated below the images) with the 
tunneling current / = 0.05 nA. 




x [nm] 

FIG. 5: (Color online) Cross sections along the longest chain 
indictated by arrow with label C in Fig. [3 for different bias 
voltages: U = -1.0, -0.75, -0.5, -0.25, 0.25, 0.5, 0.75, 1.0 
from top to bottom. The curves are shifted for better presen- 
tation. 



chains within terrace, unless we assume that gold atoms 
occupy every second site in direction a^in], which would 
be supported by fact that those chains show oscillations 
with period twice as large as the lattice constant in this 
direction. However it is difficult to find any particular 
reason supporting this assumption. So one can assume 
that we deal with one gold and one silicon chains. This 
would be supported by above-mentioned differences be- 
tween chains seen in STM images. However there might 
be another scenario in which both structures come from 
silicon whilst gold atoms substitute into the surface pro- 
ducing a third chain not visible in STM images. Such 
supposition is supported by the fact that Au is more 
electronegative than Si and binds one Si electron with 
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FIG. 6: (Color online) The length (distance between endmost 
maxima in Fig. |3| of the longest chain C indicated by arrow 
in Fig. |5|as a function of the bias voltage. 

their s, p electrons in a low laying state according 
to the first principle calculations H,0. The chain struc- 
tures visible in STM are likely to be Si atoms at step edge 
(chain C), similar to those in Si(335)-Au, and Si adatoms 
(chain D). 

III. THEORETICAL DESCRIPTION 

The experimental studies suggest that the electronic 
states in the chain C have extended nature while in chain 
D they are localized. In the language of the tight bind- 
ing model the hopping integral between atoms in chain 
C has larger value than that in the other chain or single 
particle electron energies of the chain C atoms lay much 
closer to the Fermi energy of the surface than those of the 
chain D. So one can conclude that chain C has metallic 
character whilst D is of isolating (semiconducting) na- 
ture. Moreover, the fact that at negative sample bias U 
(Fig. E|a)) chain C is more visible in the STM topog- 
raphy image than the chain D while at positive U (Fig. 
EJb)) the effect is opposite, suggests that we deal with 
mainly occupied states in the chain C and empty states 
in D. 

In order to explain the effect associated with the elec- 
tron density oscillations and differences between neigh- 
boring chains observed in STM images, we propose a 
model of two coupled non-equivalent chains composed 
of Au (chain C) and Si (chain D), schematically depicted 
in Fig. □ 

The Hamiltonian describing system composed of sur- 
face, two chains and STM tip is: 

H = HsTM + #tip + ^wire + #surf + #int (1) 

where 

HSTM = ^2 6kC k Ck ( 2 ) 

keSTM 

is for STM electrode electrons with single particle ener- 
gies 6k- The STM tip is modeled by single atom with 



FIG. 7: (Color online) Schematic view of the model STM sys- 
tem containing two chains a — 1, 2, composed of atoms with 
single particle energies £iq, intra chain hopping t XCK and inter 
chain hopping t y coupled to the surface via Vki«- STM tip is 
modeled by single atom with energy so coupled to the STM 
electrode via Vko- Parameter Uo is responsible for tunneling 
between STM tip and the chains. 

energy level 

H t [ P = £oCqC (3) 
The wire part of the Hamiltonian is described by 

Hwire — £iaCia C ia ^] t xa cf a Cj a + ^y c ta C i(^ (^) 

ia xa ia^(3 

where Si a is the atomic energy of the i-th atom in chain 
a = 1,2, t xa is the hopping integral between neighboring 
atoms in the same chain, while t y - hopping between 
neighboring atoms in different chains. Hamiltonian of 
the surface is 

#surf = ^2 6kC k C k ( 5 ) 

kGsurf 

and interactions between different subsystems are in the 
form 

Him = ^2 (^ k ° C k °0 + h' c + ^2(U0acl a C + h.C.) 
keSTM ia 

)(6) 

kGsurf ia 

with Vko being hybridization between STM electrode 
electrons and STM tip, Uo a - hopping between atoms 
in a chain a and STM tip and finally Vkia - hybridiza- 
tion matrix element connecting surface and the chain a. 
In all above equations we have omitted spin index in the 
electron creation (annihilation) operators as we are not 
interested in magnetic properties and assume that system 
is in paramagnetic phase. 

In order to calculate the tunneling current from the 
STM electrode to the surface we follow the standard 
derivations J^, H3, and get 

I= i J „ ^ T H[/STMM - /eurfM] (7) 
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where / STM ( SU rf)M = (exp ((u - /i S TM (surf)) A? 1 ) + 
is the Fermi distribution function of the STM (surf) 
electrode with the chemical potential fisTM (surf) an d the 
transmittance T(u) is given in the form 



(8) 



TstmM = 27rX)kGSTM l^ko| 2 ^(^ - ek) and r surf (cj) = 
2tt SkGsurf l^ka| 2 ^(^ — ^k) is the coupling parameter be- 
tween STM electrode and the tip atom and surface and 
the chain respectively. G r ida (uo) is the the Fourier trans- 
form of the retarded GF G r ida (t) = i0(t)([di(i), d + (0)]+), 
and is the matrix element (connecting the tip atom 
with i-th atom in chain a) of full GF, solution of the 
equation 



(ul-H)G r (u) = l 



(9) 



Full GF G r (u) is (27V + 1) x (27V + 1) matrix (TV atoms in 
a wire 1, N atoms in a wire 2 and the tip site), which has 
to be determined, just by inverting matrix (ujl — H -1 ). 
Note that within the present model there are no Coulomb 
interactions between electrons and the problem can be 
solved exactly. 



IV. NUMERICAL RESULTS 

So far we are able to calculate the tunneling current for 
a given bias voltage eU = ^stm — Msurf at constant STM 
tip - surface distance. In our model, the quantity which is 
related to this distance is the coupling between tip atom 
and the chain a Tio a = \Uo a \ 2 (in units of T = Tstm + 
r sU rf = 1)- From the physics of tunneling phenomena 
we know that Tio a exponentially depends on distance 
z between STM tip and the chain. A rough estimate of 
Tio a may be deduced by taking both the tip and the chain 
orbitals as the lowest states of square potential wells of 
spatial widths Im separated by the barrier of height equal 
to the work function and width equal to the distance z 
between STM tip and the chain. The half of the ground 
state energy splitting in this double well structure I^oa 
is thus given by |27| : 



^iOa — 



h 3 7r 2 y/2mWf 



m 2 l 3 M 



(10) 



where m is the electron mass, Wf - work function taken as 
4.5T, Im - length parameter (of order of an orbital spatial 
size) chosen as 0.52 in units of the lattice constant a^jo]- 
In numerical calculations we have chosen the Fermi 
energy of the surface as zero of energy scale (Ep — 0). 
All energies are measured with respect to it in units of 
T = Tstm + r sur f = 1 and for simplicity we have chosen 
Tstm = r sur f = 0.5. The quantities concerning distances 
are expressed in units of the lattice constant along chain 



( a [iio] = !)• Moreover, the distance between chains is 
modeled via inter-chain hooping integral t y . 

To get the information on topography of the wire we 
have solved Eqs.(0)-(0 self-consistently for a given bias 
voltage eU and fixed current /, got value of Tio a and 
finally determined z from relation (|1U|). 

In Fig. [HI we show the comparison of the STM tip - 
surface distance z along chains 1 and 2 for negative and 
positive sample biases. It is clearly seen that chain 1 
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FIG. 8: (Color online) STM tip - surface distance along chain 
1 (top panel) and chain 2 (bottom panel) for eU = —2.0 
(circles) and eU = 2.0 (squares). The other parameters are: 
en 2.0, Ei2 -0.5, e 0.0, t x 1.0 and t y = 5 • 10" 2 . 
Note that point represent positions of the STM tip (exactly 
above chain atoms) and the line is a guide for the eye. 



(corresponding to the chain C observed in experiment) 
shows the bias dependent reverse of the topography (top 
panel) while chain 2 (chain D) is insensitive to it (bot- 
tom panel). Furthermore, if the atomic energies of the 
chain are close to the Fermi energy, the topography of 
the chain remains almost unchanged with respect to the 
reverse of the bias voltage, while it reverses if the chain 
atomic energies lay far away from the Fermi energy. We 
have performed a number of calculations for various val- 
ues of Sia and this seems to be a general tendency, pro- 
vided the other parameters are the same for both chains. 
Actually this can be understood from a simple argument 
if we consider not a chain but single atom with atomic 
energy £ at . Namely, the atomic local density of states 
(LDOS) has Lorentzian shape due to the coupling to the 
STM tip and the surface, and if £ at is close to the Fermi 
energy, the LDOS is almost symmetric with respect to 
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the Fermi energy. Tunneling current is an integral of the 
LDOS bounded by the chemical potentials of the STM 
tip and surface electrodes, so it gives approximately the 
same values of the current for both voltage polarizations 
and therefore no topography reverse. In opposite situa- 
tion, if the atomic energy is far away from the the Fermi 
energy, the LDOS is strongly asymmetric and centered 
around £ a t, thus giving a different contributions to the 
current. If e at lays within eU (between /xstm and /i sur f) 
tunneling occurs through the main peak of the LDOS 
giving large value of the current while in opposite situa- 
tion the contribution is small as the tunneling takes place 
only through the tail of the main peak in the LDOS. 

Similar argument applies to the chain, however density 
of states (DOS) is now more complicated as it contains 
many peaks. Again, if the chain single particle energies 
lay close to the Fermi energy, the DOS is almost sym- 
metric and thus no topography reverse. If Ei a are much 
larger than the Fermi energy, density of states contains a 
series of peaks distributed mainly around original values 
of Si a , and therefore we get a large current for given bias 
sample and small one for reversed bias. This is clearly 
seen in Fig. where the example of the weighted DOS 
(transmittance) at STM tip position in the middle of the 
chain x = 8 (see Fig. (HJ) is plotted as a function of energy. 



is plotted for different values of the inter-chain couplings 
t y . It is clearly seen that the strength of the coupling 



1.0 



0.6 



0.4 



0.2 



o.o 



8 M =2.0 

£ i2 =-0.5 



\ — -5 1 i 1 - ' i ^ tf 



FIG. 9: (Color online) Energy dependence of the transmit- 
tance T(uj) at STM tip position x = 8 in chain 1 (solid line) 
and chain 2 (dashed line). Model parameters are the same as 
in Fig. |H1 

Turning to the experimental data, one can associate 
chain C with the chain 1 in Fig. [S| (top panel) and chain 
D with that indicated by 2 in Fig. [S| (bottom panel) and 
conclude that chain C (Fig. [2J) consists of single atoms 
with atomic energies much larger than the Fermi energy, 
while chain D has those energies close to the Fermi energy 
assuming that the other parameters characterizing chains 
are equal for both chains. 

At this point it is interesting to note boundary effects 
clearly seen in Fig. [S] (bottom panel) and associated with 
the zero dimensional states at the ends of the chains, 
recently observed on Si(553)-Au surface [25^ . 

Another important quantity is the inter-chain hopping 
integral t y . In Fig. [HI the STM tip - surface distance 
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FIG. 10: (Color online) STM tip - surface distance along 
chain 1 for eU — —2.0 a) and eU — 2.0 b) and chain 2 for 
eU = —2.0 c) and eU = 2.0 d) for different values of the 
inter- chain hopping integral t y indicated in the figure. Other 
parameters are so = 0.0, t x i = t X 2 = 1.0, en — 2.0 and 
s t2 -0.5. 

can modify results. However, while it has strong influ- 
ence on the properties of the chain 1, eventually leading 
to the reverse of the topography for positive bias (com- 
pare curve with boxes (t y = 0.8) with the others in the 
panel b)), it very weakly modifies topographic properties 
of the chain 2. Such a difference in the influence can be 
explained in the following way. As we would expect the 
inter-chain hopping t y modifies the transmittance T(ui) 
of the system. It shifts the positions of the resonances in 
the T(uj) in both chains in the same way but at the same 
time it modifies their spectral weights in different way. 
However for Si a close to the Fermi energy the transmit- 
tance integrated over the eU remains almost unchanged, 
thus leading to similar values of the tip - surface distance 
for different inter-chain couplings t y . 

Let us now discuss what the effect of the hopping in- 
tegral between between neighboring atoms in the same 
chain on the topography is. Figure ITT1 shows the topog- 
raphy of the chains for a number of the intra chain hop- 
ping integrals t X 2 with fixed t x i = 1.0 and en = 2.0, 
Si2 = — 0.5. It is clearly seen that the intra chain hop- 
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FIG. 11: (Color online) STM tip - surface distance along 
chain 1 for e£7 = —2.0 a) and eU — 2.0 b) and chain 2 for 
e£7 = —2.0 c) and eU — 2.0 d) for different values of the 
hopping integral £ X 2 indicated in the figure and fixed t x i = 1.0. 
Other parameters are eo = 0.0, t y = 0.05, en = 2.0 and 
s i2 -0.5. 



ping can change the topography of the chain. While for 
larger values of the hopping along chain 2 than that along 
chain 1 (t X 2 > t x i) one can observe reverse of the topog- 
raphy in both chains (compare curves with circles in all 
panels), for equal hoppings (t X 2 — t x i) there is a reverse 
of the topography in chain 1 and no such effect in chain 2 
(curves with boxes) and for (t X 2 < t x i) there is no reverse 
of the topography in either chain (curves with triangles). 
Such modifications of the chain topography due to the 
intra chain hopping integrals is not a general property of 
the system as there is no similar effect in the symmetric 
situation, namely if the chain atomic energies were equal. 
We have performed a number of calculations for en = Si2 
and got quite different behavior from that presented in 
Fig. Namely, our self-consistent calculations always 
converged to the results showing the reverse of the to- 
pography in both chains. This is a very important result 
as it leads to the conclusion that the chain atomic ener- 
gies have to be different in both chains in order to get an 
agreement with the experimental data. 

Similar effect on the topography has the coupling to 
the surface V\^i a . For en = e{2 one always gets reverse of 
the topography in both chains. This can be easily under- 



stood if we remind that coupling to the surface leads to 
the broadening of the resonances in the density of states 
but does not influence their positions when the Coulomb 
interactions in chains are neglected. 

From above discussion one can conclude that the main 
effect which leads to the reverse of the topography in one 
chain and its lack in the other one comes from the chain 
atomic energies - they have to be different in both chains. 
Going further one can conclude that the chains observed 
in STM have been made of different materials, one of 
silicon and the other one of gold. 



V. COMPARISON WITH EXPERIMENT 

To make quantitative comparison of the theoretical cal- 
culations with the experimental data we use the formula 
for the tip - surface distance ([TUj) and put r sur f = 1 eV 
and a = anii] = 0.384 nm. 

Figure shows both the experimental and the the- 
oretical data of typical short chains of type C and D, 
marked with smaller box in Fig. ^ The best fit has been 
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FIG. 12: (Color online) Comparison of the experimental data 
(circles) of typical short chain C for U — — 1 V a), [/ = 1 V 
b) and chain D for U = -1 V c) and U = 1 V d) with the 
theoretical calculations (squares). The model parameters are 
described in the text. 

obtained with the chain atomic energies en = 0.3 eV 
for chain C and ei2 = —0.2 eV for chain D, respectively. 
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Moreover, the values of the intra-chain and inter-chain 
hoppings were the same as those obtained in Ref. i.e. 
£xi = ^x2 = 0.68 eV and t y = 0.05 eV. Those values of the 
chain atomic energies, as it was stated previously, suggest 
that we deal with different chains, i.e. chains composed 
of different materials. Going further, we can assign the 
C chain to the Si atoms, while D chain to the Au atoms. 
This can be deduced from different values of the atomic 
energies e a used in calculations and electronegativity of 
free atoms. The Au atoms are more electronegative, so 
they should have lower atomic energies, e at = Si2 = — 0.2 
eV in our case, while Si atoms have e at = £n = 0.3 eV. 
Of course, this conclusion is true provided the other pa- 
rameters are the same or similar in both chains. 

Although theoretical topography curves do not ex- 
actly follow experimental ones, main features, i.e. the 
reverse of the topography for chain C and its lack for 
chain D, are quite reasonably reproduced. The devia- 
tions between theoretical calculations and experimental 
data steam from the fact that we used a simple model of 
STM tunneling, nevertheless giving a semi-quantitative 
agreement with the experiment. Moreover, our model al- 
lows for calculations of the tunneling current in the cases 
when the STM tip is placed exactly above the chain atom, 
therefore we do not present the theoretical data for the 
tip position between chain atoms. To get full agreement 
with experimental data one has to take into account ad- 
ditional effects, like tunneling into aside chain atoms, the 
Coulomb interactions or perform first principles calcula- 
tions which is out of the scope of the present work. 



VI. CONCLUSIONS 

In conclusion we have studied properties of the mono- 
atomic chains on the vicinal Si(557)-Au surface using 
STM technique. We have observed double non-equivalent 
chains on such surface showing different behavior of the 
topography with reversing of the bias voltage. While 
chain C shows reverse of the topography, chain D does 
not show such effect. Additionally the STM tip - sur- 
face distance oscillates along any chain with the same 
period equal to double lattice constant in direction [110]. 
We have also performed tight binding model calculations 
showing similar oscillations of the STM tip - surface dis- 
tance and led us to the conclusion that atomic energies 
have to be different in different chains and larger in chain 
C than in chain D. This fact supports the scenario which 
is probably realized in experiment that one chain is made 
of gold (chain D) and the other one of silicon atoms (chain 
C). 



Acknowledgments 

This work has been supported by the KBN grant no. 
1 P03B 004 28. T. K. thanks the Foundation for Polish 
Science for a Fellowship for Young Scientists. 



[1] F. J. Himpsel, K. N. Altmann, R. Bennewitz, J. N. Crain, 
A. Kirakosian, J. -L. Lin, J. L. McChesney, J. Phys.: 
Condens. Matter 13, 11097 (2001). 

[2] J. M. Luttinger, J. Math. Phys. 4, 1154 (1963). 

[3] F. D. M. Haldane, J. Phys. C: Solid State Phys. 14, 2585 
(1981). 

[4] R. E. Peierls, in Quantum Theory of Solids, Clarendon, 
Oxford (1955). 

[5] M. Jalochowski, E. Bauer, Prog. Surf. Sci. 67, 79 (2001); 

Surf. Sci. 480, 109 (2001). 
[6] R. Zdyb, M. Strozak, M. Jalochowski, Vacuum 63, 107 

(2001). 

[7] J. N. Crain, A. Kirakosian, K. N. Altmann, C. 
Bromberger, S. C. Erwin, J. L. McChesney, J. -L. Lin, F. 
J. Himpsel, Phys. Rev. Lett. 90, 176805 (2003). 

[8] J. N. Crain, J. L. McChesney, Fan Zheng, M. C. Gal- 
lagher, P. C. Snijders, M. Bissen, C. Gundelach, S. C. 
Erwin, F. J. Himpsel, Phys. Rev. B69, 125401 (2004). 

[9] M. Krawiec, T. Kwapinski, M. Jalochowski, phys. stat. 
sol. (b) 242, 332 (2005). 
[10] M. Henzler, R. Zhachuk, Thin Solid Films 428, 129 
(2003). 

[11] M. Jalochowski, M. Strozak, R. Zdyb, Surf. Sci. 375, 203 
(1997). 

[12] P. Segovia, D. Purdie, M. Hengsberger, Y. Baer, Nature 

(London) 402, 504 (1999). 
[13] R. Losio, K. N. Altmann, A. Kirakosian, J. -L. Lin, D. 



Y. Petrovykh, F. J. Himpsel, Phys. Rev. Lett. 86, 4632 
(2001). 

[14] K. N. Altmann, J. N. Crain, A. Kirakosian, J. -L. Lin, 
D. Y. Petrovykh, F. J. Himpsel, Phys. Rev. B64, 035406 
(2001). 

[15] D. Sanchez-Portal, J. D. Gale, A. Garcia, R. M. Martin, 

Phys. Rev. B65, 081401 (2002). 
[16] I. K. Robinson, P. A. Bennett, F. J. Himpsel, Phys. Rev. 

Lett. 88, 096104 (2002). 
[17] D. Sanchez-Portal, R. M. Martin, Surf. Sci. 532, 655 

(2003) . 

[18] H. Okino, R. Hobara, I. Matsuda, T. Kanagawa, S. 

Hasegawa, J. Okabayashi, S. Toyoda, M. Oshima, K. 

Ono, Phys. Rev. B70, 113404 (2004). 
[19] R. A. Zhachuk, S. A. Teys, A. E. Dolbak, B. Z. Olshanet- 

sky, Surf. Sci. 565, 37 (2004). 
[20] A. A. Baski, K. M. Jones, K. M. Saoud, Ultramicroscopy 

86, 23 (2001); A. A. Baski, K. M. Saoud, K. M. Jones, 

Appl. Surf. Sci. 182, 216 (2001); J. W. Dickinson, J. C. 

Moore, A. A. Baski, Surf. Sci. 561, 193 (2004). 
[21] S. S. Lee, N. D. Kim, C. G. Hwang, H. J. Song, J. W. 

Chung, Phys. Rev. B66, 115317 (2002). 
[22] J. R. Ahn, Y. J. Kim, H. S. Lee, C. C. Hwang, B. S. Kim, 

H. W. Yeom, Phys. Rev. B66, 153403 (2002). 
[23] S. Jeong, H. Jeong, S. Cho, J. M. Seo, Surf. Sci. 557, 183 

(2004) . 

[24] S. C. Erwin, Phys. Rev. Lett. 91, 206101 (2003). 



9 



[25] J. N. Crain, D. T. Pierce, Science 307, 703 (2005). 

[26] S. G. Davison, M. Steslicka, Basic Theory of Surface 

States, Oxford Univ. Press, New York (1992). 
[27] Y. Calev, H. Cohen, G. Cuniberti, A. Nitzan, D. Porath, 

Israel J. Chem. 44, 133 (2004). 
[28] H. Haug, A. P. Yauho, Quantum Kinetics in Transport 

and Optics of Semiconductors, Springer, Berlin (1996). 



[29] L. V. Keldysh, Sov. Phys. JETP 20, 10 108 (1965). 
[30] C. Niu, D. L. Lin, T. -H. Lin, J. Phys.: Condens. Matter 

11, 1511 (1999). 
[31] M. Krawiec, K. I. Wysokihski, Solid State Commun. 115, 

141 (2000); Supercond. Sci. Technol. 17, 103 (2004). 



